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ABSTRACT 


A mathematical model is developed to determine the heat 
transfer rate of a cylindrical condenser section of a rotating 
heat pipe. This model is coupled to an existing code and an 
analysis is accomplished on both a smooth and axially finned 
Sondensers, |thelresults of this analysis are compared to those 
Of a Samilar analysis performed on a tapered condenser heat 
pape using adentical geometric and operating parameters. 

Results of the comparison indicate cylindrical condensers 
have a lower heat transfer rate than equivalent tapered,con- 
Semsenrs. sllis reduction ianvhéeat transfer rate is dueyto a 
greater condensate film thickness and is most significant in 
a smooth condenser. 

Axially finned condensers with triangular and rectangular 
fin profiles are also compared. The rectangular fins are 
assumed to have adiabatic tips. Results indicate the heat 
Eransfer rates for these two profiles vary by only 0.40 per 


cent for both tapered and cylindrical condenser. 
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I. INTRODUCTION 


me THE ROTATING HEAT PIPE 

ihe rotating, wickless heat pipe is a closed container 
designed to transfer large amounts of heat from rotating 
Machinery components. Essentially, it consists of three 
Maan components: a cylindrical evaporator section, a con- 
denser section which may be either tapered or cylindrical in 
shape, and a fixed amount of working fluid. A typical tapered 
rotating heat pipe is; shown in Figure 1. 

When the heat pipe is rotated about its longitudinal axis 
mea speed above a certain critical value, the working fluid 
forms an annulus in the evaporator section. Note in Figure l 
eeat the diameter of the evaporator is larger than the con- 
Memser. This larger diameter provides a, greater liquid reser- 
weer. AS heat is added’ to the evaporator, the fluid in the 
evaporator will vaporize. The vapor will flow axially towards 
the condenser as a result of a slight pressure difference, 
imansporting' the latent heat of vaporization with it. In the 
condenser end, external cooling of the condenser causes the 
WapOr to condense. In the case of a tapered heat pipe, the 
Gemtritugal force due to the rotation of the pipe has a con- 
ponent acting along the condenser wall which accelerates the 


liquid condensate back to the evaporator to complete the cycle. 
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Schematic Drawing of a Tapered 
Condenser Rotating Heat Pipe 
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Mvicylindrical heat pipe, on the other hand relies on a hydro- 
static pressure gradient to drive the liquid condensate back 


to the evaporator. 


B. BACKGROUND 

The first theoretical investigation into the performance 
of a tapered rotating heat pipe at the Naval Postgraduate 
menoo] was accomplished by Ballback [Ref. 1] in 1969. He 
examined the limits in heat transfer controlled primarily by 
mimad dynamic considerations. In particular, he considered 
|me following four limits on heat pipe performance: a) the 
Gerling limit, b) the entrainment limit, c) the sonic limit 
aad ac) the condensing limit. Tantrakul [Ref. 2] calculated 
@mese limits for a specific heat pipe. He found the condens- 
meee limit was the controlling limitation. In fact, the cal- 
culated heat transfer rate, based on the condensing limit was 
im eOth the heat transfer rate for the next: lowest limit, the 
entrainment limit. 

in order to overcome this.condensing limitation and thus 
mmerease the heat transfer rate of the rotating heat pipe, 
Emme concept of an internally finned tapered rotating heat 
pepe was considered by Schafer [Ref. 3]. Schafer developed 
a@m@eanalytical model for this tapered heat pipe with a trian- 
gular fin profile as shown in Figure 2. He assumed one dimen- 
stondal heat conduction through the wall and fin. Corley [Ref. 


4] developed a two-dimensional heat conduction model using a 
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Figure 2. Axially Finned Condenser Geometry Showing 
Rims, .noughs, and, Lines. of Symmetry 


18 


Finite Element Method formulation for the same geometry. To 
overcome the problem of few nodal points along the fin surface, 
Corley assumed a parabolic temperature profile along the fin 
Surtace. Tantrakul [Ref. 2] modified Corley's computer code by 
imereasing the number»sof finite elements from two to three in 
@eder to minimize the error at the apex of the fin. Purnomo 
[Ref. 5] developed a two-dimensional Finite Element Method 
solution to the steady state heat conduction problem using a 
Minear triangular finite element. Davis [Ref. 6] modified 
Purnomo's code to make it compatible with COPES/CONMIN [Ref. 8], 
aoe Optimization program. Davis, in his modification, found a 
Seding error in Purnomo's [Ref. 5] code, that once corrected, 
Permitted Purnomo's corrected code to converge to Schafer's 
[Ret. 3] results. 

Purnomo’s {Ref. 5] code is limited in that it is restricted 
memone particular condenser geometric configuration, namely: 
an axially finned tapered condenser heat pipe with a triangular 
oem profile. In that tapered finned condenser are difficult 
to manufacture, it is doubtful that any widespread practical 
aeplication of this geometric configuration will result. Cylin- 
drical condensers, on the other hand, can be manufactured with 
meen tess difficulty and might find practical application. 

This being the case, a more practical and beneficial code 
would be one that could analyze cylindrical condenser rotating 
Neat pipes, both finned and smooth. In actuality, the most 


beneficial code would be one that could analyze the following 
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mou geometric CONtigurations: 1) tapered-internally finned, 
2) tapered-smooth, 3) cylindrical-internally finned, and 4) 
Gylindrical-smooth. The theoretical heat transfer performance 
of the four geometries could then be compared to determine the 
advantages and disadvantages of each design. An additional 
advantage would be gained if different fin profiles, i.e., 
triangular vice rectangular, could also be analyzed and 


compared. 


See THESIS OBJECTIVES 
The objectives of this thesis are: 

1) Develop analytical models for both cylindrical-smooth 
and cylindrical-axially finned condensers. 

2) Develop solution techniques to these analytical models 
that will account for temperature variations along the 
axial length of the condenser. 

—eeModiry Purnomo'’s [Ref. 5] code to provide a solution to 
the two-dimensional steady state conduction heat transfer 
problem for the following four geometric configurations: 
a) tapered-smooth, b) tapered-finned, c) cylindrical- 
smooth, and d) cylindrical-finned. 

4) Modify Purnomo's [Ref. 5] code to provide the additional 
capability Of analyzing a rectangular fin profile with 
an adiabatic tip. 

5) Obtain and compare results of the four geometric con- 


figurations given above for various operating conditions. 
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Ti. THEORETICAL ANALYSIS FOR A CYLINDRICAL HEAT PIPE 


A. INTRODUCTION 

in at cylundrical condenser heat’ pape, "the radius of the 
condenser is constant along the axial length of the condenser. 
The flow of the condensate in the absence of vapor-liquid 
ioterfacial shear, is dependent upon the variation in hydro- 
Static pressure with changes in film thickness along the sur- 
face of the heat pipe. Leppert and Nimmo [Refs. 8 and 9] 
investigated the phenomenon of film condensation on a flat 
Merazontal plate. This situation is similar to film conden- 
meron On the inside surface of a rotating cylindrical con- 
demser. In the case of a cylindrical condenser, the body 
mimee, rather than being the force of gravity, is now the 
Sentrifugal force caused by the rotation of the heat pipe. 
Weigenseil [Ref. 10] and Tantrakul [Ref. 2] compared experi- 
memeat results for a cylindrical condenser rotating heat pipe 
with the theoretical results of Leppert and Nimmo [Refs. 8 
and 9] and found good agreement. The Leppert and Nimmo solu- 
tion was limited in that it was based on a constant surface 
Ecmperature along the length of the plate. A rotating heat 
pipe, in actuality, has a temperature variation along the 
axial length of the condenser which in some cases, may be 
Significant. This being the case, it was necessary to develop 


a mathematical model which would consider the axial temperature 
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faeracion ingthe.solution of the heat transfer :analysis., In 
the mathematical development that follows, a cylindrical 

smooth (unfinned) condenser will first be considered in that 
it is the simplest case. The model will then be extended to 


include a cylindrical axially finned condenser. 


B. THEORY FOR A CYLINDRICAL SMOOTH CONDENSER 
1. Assumptions 
In developing the theoretical analysis, the following 
assumptions are made: 
a) Film condensation, not dropwise condensation occurs in 
the condenser. 
b) The condensate film undergoes laminar flow. 
c) Momentum changes through the condensate are small. 
Thus, there is essentially a static balance of forces. 
aye ihe vapor exerts no drag in the condensate; there is 
no interfacial shear. 
e) The temperature distribution within the film is linear. 
f) The vapor space is essentially at one pressure, gas 
g) The density of the fluid is much greater than the den- 
Saity Of the-vapor. Thus, the density*of the vapor can 
be neglected. 
fh) The centrifugal force is much greater than the force 
of gravity and, thus, gravity may be neglected. 
i) Velocity gradients in the circumferential direction 


relative to the pipe wall are negligible. 
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ie the condensate film thickness is much less than the 
radius of curvature of the condenser wall. 
k) The rotating heat pipe is operating at steady state 
conditions. 
2. Condensate Momentum Equation (X-Direction) 

By applying the above assumptions and the coordinate 
system shown in Figure 3, an analysis similar to Nusselt's 
eraginal film condensation film theory may be used [Ref. 11]. 
Based on assumption c, a static force balance may be taken 
Gmaean infinitesimal fluid element in the x-direction as shown 
iaerieure 3. This force balance results’ in the following 


equation: 


ZF Y =i Qe Ley Jame 0 (eqn 2.1) 


where t = shear stress (1bf/ft? ) 
p = pressure (lbf/ft? ) 
= co-ordinate Measuring distance aliong surface (ft). 
y = co-ordinate measuring distance normal to surface (ft). 
3. Condensate Momentum Equation (Y-Direction) 
inva samitar manner, Using Figure 3; a force balance 
mm the y-direction yields: 
Be .= 0 : a P_w'r = 0 (eqn, 2.2) 


where I density of the fluid (lbm/ft?) 


w angulbandvelocitys(rad/hr ) 
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Figure 3. Cross Section of Infinitesimal Fluid Element on 
Cylindrical Condenser Intetnal Surface 
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Integrating equation (2.2) between the limits y and 
é* for y and corresponding limits of P and P, for pressure 


results in the following equation: 


Pp = Bees OP pw°r (o*=¥) (eqn 2.3) 


where ue = the pressure of the vapor (lbf/ft*), and 


§* 


film thickness (ft) 
Differentiating equation (2.3) with respect to x yields the 


following expression for dP/dx: 


Hee dPy 5 eee «66 
dx dx dx 


(eqn 2.4) 


applying assumption (f), (P, is constant, therefore, dP, /dx=0) 


ei@e substituting equation (Z.4) into equation (2.1) yields: 


ay = p -w°r ae (eqn Z..5) 


Integrating equation (2.5) with the corresponding limits of 
Ba@eecration y to 6* and t to 0 results in the following ex- 


pression for shear stress: 


pe dé* 
dx 


ye or | (eqn 2.6) 
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eS U ou (eqn 2.7) 


where UW = fluid dynamic viscosity (lbm/ft-hr) 


condensate velocity (ft/hr) 


u 
Substituting equation (2.7) into equation (2.6) and integrating 


with the corresponding limits of integration 0 to y and 0 to u 


yields: 
2 
Ow Tr = 2 
d c 
u = ——— = [ - - yé*] (eqn 2.8) 


The average velocity of the condensate may be found in the 


following manner: 


* — 
b : ) 1 p pwr dé* 2 : 
5 = Ei udy = 4 f [ae [ay (eqn 2.9) 
6 fe) 6 fe) u ax = 


or 
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i) 
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(eqn 2.10) 


a Continuity Equation 


the CcoOmtinuityrequation for mass. flow’ requires® that: 


h = P¢ uA (eqn 2.11) 


where m™m condensate mass flow rate (lbm/hr) 


Mo Mcross seetronal area of the fluid (£t2) 
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This can. also be written as 


* 


6 ahs 
m Ih Pp .u 2mrdy (eqn 2.12) 
O 


aubstituting equation (2.10) into equation (2.12) and inte- 


grating yields: 


72 
2,2 * * 3 
m - .2T0f wir - a (eqn 2418) 
u Xx 


Differentiating this equation with respect to x yields: 


do, d6* 5%" 
dx u dx dx 3 


] (eqn 2.14) 


6. Energy Equation 


Having applied assumption (e), if the film surface 
temperature is at the saturation temperature (tae) of the 
Wapor and if the wall of the axial increment is at a given 
constant temperature (Ts then the heat transfer by conduc- 


con of a fluid element of surface area dA is: 


Ke (Tsat - Tw) dA 
a ee 


dq = : (eqn 2.15) 
where dq = differential heat transfer rate (Btu/hr) 

dA = 2nrdx(ft7) 

ky = thermal conductivity of the condensate film 


(Btu/ar=£t-=F ) 
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la 
i) 


nee saturation temperature (degrees F) 


J 
i] 


inside condenser wall temperature (degrees F) 

Considering the change of phase and defining Ags as 
the average enthalpy change of the vapor in condensing to a 
liquid and subcooling to the average liquid temperature of 


the film, then dq is also defined by: 


dq = Agy di (eqn. 2..10) 
where he, = laten heat of vaporization (Btu/1lbm) 
Cy = specific heat (Btu/1bm R) 
AT = (Tsar - Ty? 
eae hi 


fg $2°05-55.—2 Cy Sa AT, 


Rearranging equation (2.16) and substituting this equation into 


equation (2.15) yields: 


oe 6 (éqn-: 2.17) 


Panally coupling the energy and continuity equations result in 


the following differential equation: 


ry d dsé* * 3 4 a ) 
eet tin, a2 Ske (sat = Tw)u Cs Tw) u Geant 2918) 
Pew Rey 


Equation (2.18) can be solved using the Finite Element 


Method to provide the film thickness profile along the axial 
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iength of a cylindrical condenser. Appendix A provides a 
Peeatileducgescripeion o£ thisssolution. Once the film pro- 
file is known, a steady state two-dimensional heat conduc- 
tion analysis can be performed. 
7. Determination of Heat Transfer Rate 

Assume that the cylindrical condenser section of the 
mOtating heat pipe is divided axially into a number of incre- 
ments. Then for any axial increment of a cylindrical conden- 
Ser, the differential heat flux can be determined by the 


following expression. 


i = Ah ee (eqn Zk 9} 
Ora, Chick; 41 
Ke Ky Reet 
where To = ambient temperature (degrees F) 


thick = thickness of the condenser wall (ft) 
k = thermal conductivity of the wall material 
(Beu/hr- ft -.) 
DN ack = external heat transfer coefficient 
(Btu/hr-ft?-F) 
Note the three terms in the denominator are the thermal re- 
Sistances of the film, wall and external convection respectively. 


the differential heat transfer rate for any. increment 


can be found by the following relationships: 


dq = dq’! + 2nrdx (eqn 2.20) 
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(eqn 21.1271) 


Bamacion (2.21) represents the total heat transfer rate for 
maeincremental section of width dx. To find the total heat 
feanster rate for the entire cylindrical condenser, the in- 
cremental heat rates must be summed over the entire length 
Of the condenser. Therefore: 

NDIV 


2 dq (eqn 2.22) 


Qtotal re 


where NDIV = total number of axial increments. 


€. THEORY FOR A CYLINDRICAL AXIALLY FINNED CONDENSER 
i ASSUMptiONs 
Referring to Figure 4, it is obvious that the analysis 

Of a cylindrical internally finned condenser is more compli- 
cated due to the mass flow from the fins into the trough 
region between the fins. For this reason, in addition to 
the simplifying assumptions made for the smooth condenser 
match are listed in the previous section, the following 
assumptions must also be made: 

mee Referring to Figure 5, the mass flow along the fin sur- 


face does not flow axially in the x-direction, but only 


30 


12ST Get 


lo pelbisse Tg 


Hs Fa oie" 
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Cc) 


aloneethe sumtace Of the fin im the z-direction into the 
trough. Thus, mass tlow in the axial direction is only 
permitted in the trough region between the fins. This 
is a reasonable assumption in that the film thickness 
along the fin surface is very small in relation to the 
flim Badekwess in the trough. This being /the case, the 
nyadrostatic tomnec in the x-direction onfa fin fluid ele- 
ment wiil*be, much less than the centrifugal force com- 
ponent) in the \z-direction ‘on that same fluid element 
forcing that tiuid elementiinto the trough. 

VusSt as an the axial direction, there is no pressure 
eEhangeyalone ithe Surface of the fin”im the z-direction. 
ie wilt be assumed that the temperature along the con- 
ywectave Surtace Of the fin 1S at a constant value (Tayg)> 
Dies average fin Surface temperature 1s the arithmetic 
menage Of the fin tip temperature and the fin base sur- 
face temperature where the fin intersects with the wall 
of the condenser. This is a valid assumption if the fin 
Section is divided into a sufficient number of finite 
elements. Purnomo's [Ref. 5] results indicate a less 
than one degree variation, even for very large fin half 
angles. This variation in temperature will have an in- 
Significant effect on film thickness along the surface 
of the fin and can be neglected by using an average 


value. 
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mesure 4.° (Trough Section of Axially Finned Condenser 


S72 


ve — 
ae 


Figure 5. Axially Finned Condenser Symmetric Fin Section 
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2. Mass Flow in the X-Direction 
As a result of assumption (a), the resulting momentum 
@€quation in both the x and y directions as well as the equa- 
tions for velocity and mean velocity are identical to those 
developed in section B for the smooth condenser and will not 
be redeveloped here. Looking now at mass flow in the x-direc- 
mien Which 1s limited to flow only in the trough, the mass 


flow rate is given by the following expression: 


Zone * 
; eite Or. ds a ee 2 2 
m S'5 Ceo" +56" tana) (ean 2.235) 


where a £in-halfoangle t(radians) 


E Width of the trough (£€) 

Nete that the quantity in parentheses is the cross sectional 
area of the film condensate in the trough (See Figure 4). 
iaeang the derivative of equation (2:23) with respect to x 


feetds*the rate of change of mass flow in*the trough for a 


given axial increment. 


dim OP -w*r x 
pe = aac & [$2 (e6*"+ 6**tana)] (eqn 2.24) 
x u 


Equation (2.24) represents the rate of change of the 
total mass flow rate with respect to x in the x-direction. 
This equation must be coupled with the energy equations for 
the fin and trough to develop a representation of the film 


peetile in the trough. 
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Ss.) Mass Plow im the Z-Direction 
Examining dn infinitesimal fluid element on the sur- 
Hace Of a fin for any axial increment of width Ax, as shown 


in Figure 6, the momentum equation in the z-direction becomes: 


— - Pew rcosa (eqn 2.25) 


where c Sliear stress imsthe zodirection (1bf/ft7) 


N 
i) 


co-ordinate measuring distance along the surface 
Of the fin: (ft) 
Neglecting dP/dz based on assumption (b), and integrating 


Equation (2.25) from tT, tO QO and y to 6 yields: 


T = ——— p pw rcosa(6-y) Geqn 2.26) 


where w Fluid velocity. in the z-direction (ft/hr) 


[e7) 
it} 


Pin vteim thickness alone the surface of the fin (ft) 
Nece, 6 , the fin film thickness should not be confused with 
6, the film thickness in the trough. Integrating equation 
(2.26) from 0 to w and 0 to y provides the following expression 
mom fluid velocity: 


p,d*rcosa 2 
Sy = - ) (eqn 2.27) 
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Figure 6. Cross Section of Infinitesimal Fluid Element on 


Pinesurtace of Axially Finned Condenser 
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This relationship may be used to find the average fin fluid 
velocity w: 
w?rcosa 6(z)? 


6 0 
weif ee me (eqn 2.28) 
6 %0 Sul 


ime mass flow rate in the z-direction along the surface of the 


PmetOr a given axaal increment is given by: 


they = p pwdA (eqn (2.29) 
where dA = the cross sectional area of the fluid flowing along 
the fin ‘surface’ (ft* ) 
Sbstituting equation (2.28) into equation (2.29) yields: 
Pw rcosa 6(z) dx : 
alee i =e Ceqn> 2°30) 
irs equation is identical to equation (15) of Schafer's [Ref. 
3] analysis if the condenser cone half angle (f) is set equal 
Bon0 which is the case for a cylindrical condenser. 
a enerey Bquation for the Trough Condensate 
An energy balance’ on an infinitesimal fluid element in 
the trough of an axial increment of width dx with surface area 
Ewrdx yields the following expression for heat transfer by 


conduction: 


44+ rough ee (eqn 2.31) 
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Note also, that the trough heat transfer rate is given by: 


dd trough 7 Nes 2 


i, rough (eqn 2.32) 


Combining equations (2.31) and (2.32) and dividing by dx re- 


sults in the following: 


dm k,-(T ive 
trough _, £ “22. W (eqns 2.55) 
dx Reg 6 


E@m@ation (2.33) is an expression for incremental change in mass 
Blow rate with respect to x due to condensation in the trough 
region. 

5. Energy Equation for the Fin Condensate 

An energy balance on a differential element of surface 

ameed dxedz yields the following relationship for differential 
Meat into the fin: 
Re oat - T 53 (2)) dxdz 


elie ae se ee (eqn 2.34) 


OG >. — = 
£in Fin 5 (z) 


where T (z) = fin surface temperature at some position z 


tan 
along the surtace, of the. fin, (degrees _F) 

Since the fin condensate mass is assumed to flow only in the 

Zeaarection, equation (2.30) is differentiated with respect 


Eee. and, substituted into equation (2.34). After substitution 


and rearrangement, the following equation results: 
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sat * Tein (2) )dz 


6(z)*dé(z) = (eqn 2.55) 
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Applying assumption (c) i.e., Tei, (2) equals Tavg fOr 


mime z and integrating equation (2.35) from 0 to 6 and 0 to z 
weelds the following relationship for fin film thickness 5 (z): 


i/s 


4 kel eat - ee 


CZ) = | (eqn 2.36) 


ae 
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where Tavg = average fin surface temperature (degrees F). 
Suestituting equation (2.36) into equation (2.30) and solving 
Bemerate Of change of mass flow rate of the fin with respect 


Pees fOr an increment of width dx yields: 


1/4 
din _- Zion w- reosa «4, kaCT OT eee 
fine te rin t | fiesat avg (an 2.37) 
QZ: 
dx 3u Op wr ee cosa 
mmere zZ* = z-¢*/cos(a). Note, z* is the distance along the 


mimeace Of the fin from the fin tip to the trough film thick- 
ness (5*). Note also, that the right hand side of equation 
Wo?) 1S multiplied by two; this accounts for mass flow from 
Pee fins on both sides of the trough. 
6. Continuity Equation 
For any axial increment of length dx, continuity dic- 


Paces that: 
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«1 (5) noitameean? 


din Gin a6 dim 
COL ad _ fan trough 
SSS Se = = (eqn 2.38) 


Eeestituting equations (2.24), (2.33) and (2.37) into equation 


(2.538) and rearranging yields: 


a = aS feo*° 2 6* tana) Js"= - 


sat avg 


i k,(T =. yz 
— (eqn: 2.39) 


- 26*cosa 
Pp w'r ey cosa 


Equation (2.39) can be solved using the Finite Element Method 
formulation provided in Appendix A. The solution of this equa- 
fron provides the film thickness profile along the axial length 
@2ea cylindrical finned condenser. 
7. Determination of the Heat Transfer Rate 
Once the film’ profile has been determined within the 
mean, the local convective heat transfer coefficient can be 


found for the trough using the following relationship: 


Ke 
Bs) oush,. Eee (eqn 2.40) 
§ (x) 


In a similar manner, the local heat transfer coefficient along 


the surface of the fin can be found by: 


k 
af 
h(Z) ein = S(z) (eqn DEEN 
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Mire ditferential heat transfer rate for any fin section, as 
Enown in Figure 4, of axial incremental length dx is: 


ay) for the trough: 


foot Tin ode Goundak 


where €°*dx is the surface area of the trough, and 


Poe tor the fin surface: 


Zo (T Pee xy SGC Z 

dd-; = “adh pia aes * (eqn 2.43) 

in 
e A(Z)e¢in 

where zis the surface length of the fin. 

ieemeotal differential heat transfer rate per axial increment 

ws tound by summing equation (2.42) and (2.43) for the total 


mumper Of fins. That is: 


NFIN 
dq otal - ; Geta : dq+ rough? (eqn 2.44) 


where NFIN is the total number of axial fins. 
iaea Similar manner, the total heat transfer rate for the entire 


finned condenser can be found by the following relationship: 


NDIV 
Qtotal . ‘ dq otal (eqn 2.45) 


where NDIV is the total number of axial increments. 
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A. 


Tit, COMPUTER. CODE DESCRIPTION 


GENERAL DESCRIPTION OF CODE 


The computer code consists of a main body and eight sub- 


routines. Basically, the code which is provided in Appendix 


fers a modification of Purnomo's [Ref. 5] code. The function 


of each subroutine used in the code is as follows: 


a) 


b) 


€) 


d) 


"CORRES" established the correspondence between the 
local and global nodal points used in the Finite Element 
Method solution for the two-dimensional steady state 
heat conduction-problem. In so doing, "CORRES" also 
numbers all elements and nodal points inthe finite ele- 
ment model and assigns local nodal points to each of the 
elements. In -additiaon "CORRES" also defines. major ele- 
ment numbers used in other subroutines as control 
parameters. 

SCOORD” defines the x and y coordinates for all- nodal 
points in the finite element heat conduction problem 
model. 

PDLSTAR" determines the film thickness (§*) on the sur- 
face of a smooth condenser or in the trough in a finned 
condenser. 

SHICOBr’ determines the heat transfer coefficient for 


all convective surface elements. 
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e) "FORMAF'" formulates the Finite Element Method equations 


for the two-dimensional steady state heat conduction 


problem. 
f) '"BANDEC" is an equation solver for a symmetric matrix 
which has been transformed into banded form. ''BANDEC" 


will return the solution to the two-dimensional heat 
conduction problem. 

eyo)" HTCALC™ determines the elemental, incremental. and. total 
heatetranster rates. 

my) ~"DELCRV" determines the condensate film profile in a cyl- 
indrical. condenser: 

Two additional Naval Postgraduate School computer library 
Toutines are also used in the code: 

ime DPOLRIY:. is a.nonIMSL. double precision library routine 
that determines the roots of a real polynomial. This 
routine is called by "DLSTAR" to determine the film 
thickness for the succeeding increment in the analysis 
Gi a tapered condenser. 

fy LEQTZF’ is an IMSL double precision libary routine that 
solves a set of simultaneous linear equations. This 
ZwOutine 1s called by "DELCRV" to solve the Finite Element 
Method equations for the cylindrical condenser film pro- 
Ele problem. The resulting film profile is: then used 
ticthe heat conduction analysis. 

im order to use the computer code to analyze heat transfer 


in a rotating heat pipe, nine data cards are required. A user's 
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guide describing these data cards and required input is pro- 
vided in Appendix B. The input data, describing the geometric 
eenfiguration of the rotating heat pipe as well as the operating 
parameters determines which solution technique is utilized in the 
meelysisy o)Thersolutionstechniquesfor each of the four con- 
denser geometries, i.e., tapered-smooth, tapered-axially finned, 
eyeindricai-smooth and icylindrical-axially finned is different. 
In all cases however, the Finite Element Method is used to 

solve the two-dimensional steady state heat conduction problem. 
This solution is the one developed by Purnomo [Ref. 5] and 

has not been modified. Details of the development of this 
solution are described in detail in Purnomo's thesis [Ref. 5] 
euamwillinot be repeated:here. This being the case, each of 


the four solution techniques will now be discussed in detail. 


B. INTERNALLY FINNED TAPERED CONDENSER SOLUTION 

The complete development of this solution technique is 
provided in Purnomo's [Ref. 5] thesis and will not be rede- 
veloped. When an equation is required for clarity, the equa- 
tion in final form will be provided. Where there is a 
modification to Purnomo's [Ref. 5] code, this modification 
will be noted. 

The condenser of the internally finned-tapered condenser 
is divided into NDIV axial increments. These axial incre- 
ments are then subdivided circumferentially into ZFIN number 


of subincrements where ZFIN is the total number of fins. 
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These subincrements are then divided in half to form the 
Basic symmetric unit, one-half of a fin-trough section as 
Suown in Figure 1. This unit section is then divided into 
a number of linear triangular elements. The number of ele- 
ments depend on the input parameters. The only limitation 
is that the number of system nodal points must not exceed 
mOOQ. otherwise, certain variables, e.g., x and y, would 
exceed allotted storage values. Figure 7 shows a unit 
Section subdivided into 25 elements. After this unit is 
Subdivided, each nodal point is assigned an x and y coordi- 
nated based on the geometric input parameters. 

Mo Start jthe/ iteration; two initial values are required: 
1) an initial temperature for the nodal points along the 
amternal convective boundary, and 2) an initial trough film 
thickness (6*). The initial temperature is provided as an 
input parameter and the initial trough film thickness at the 
feost increment.is provided by a relationship taken from an 
analysis by Sparrow and Gregg for condensation on a rotating 
disk [Ref. 12]. 

Once these values are known, the heat transfer coeffi- 
cient for the internal convective surface elements are found 


using the following relationships: 


ke Keo ¢°w* €+xsinf he cosPcosa 


ae... - = | Ceqn: 3.1) 
eam = § (Zz) 4u(-AA+z /3-BB+z"/2+(T,,,-T)z 
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Fisure 7. Axially Finned Condenser Symmetric Section Sub- 
divided into 25 Binear Triangular Finite Elements 


distance from condenser end (x=0.0) to midpoint 


miaere x 
of increment (ft) 
ta = fin apex temperature (degrees F) 
AA, BB = constants in the parabolic temperature determined 
by a Langrangian fit. 
mer the trough’ surface elements: 
Ke 
A(X) + rough = ee) (eqn 3.2) 

These heat transfer coefficients, along with the thermal 
conductivity and x and y coordinates of the nodal points are 
used to form the Finite Element Method equations for the two- 
dimensional heat conduction problem. The equations are then 
Seived to yield a temperature distribution in the symmetric 
Section. 

The above iteration is repeated, where, now the solution 
Eemperature distribution from the previous iteration is used 
memcaiculate the heat transfer coefficients along the con- 
vective surface fin elements as well as a new 6* using Sparrow's 
maeGrege's relationship [Ref. 12]. This new 6*, in turn, is 
moea tO determine the heat transfer coefficients of the trough 
elements. 

Again, the Finite Element Method subroutines will yield 
mecemperature distribution for the symmetric section. At 
this point, the nodal point temperatures are checked for con- 


vergence using the following relationship: 
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where NSNP number of system nodal points 
j = present iteration 


j-l = previous iteration 


Cie convergence criterion 

Prethis convercence test.1s successful, for,ali.nodal points, 
the increment is considered solved. If any one nodal point 
mes. the iterative process is.repeated until convergence is 
met. The convergence test in equation (3.3) is different than 
the one used by Purnomo [Ref. 5] in his thesis. Purnomo com- 
fmemed incremental heat transfer rates per. unit, of condenser 
length, Qae rather than temperature as is done in the modified 
code. 

If convergence is met, the heat transfer rate is determined. 
From this heat transfer rate, the incremental mass flow rate is 
determined by the following equation: 


2Q;Ax 


Oral. Piva (eqn 3.4) 
oO 


where Q; heat transfer rate per unit length (Btu/hr-ft) 
Ax = incremental width (ft) 
Using this value of incremental mass flow rate determined 
Pyeequation (3.4), the following equation is used to calculate 


Bie subsequent interval's trough condensate film thickness (¢*) 


With a polynomial rootfinder subroutine: 
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A pw (r+xsinf) 6** (x)sing(6* (x) e+6** (x) tana) 


m = 
total 3u 


(eqn 3.5) 


where € = trough width (ft) 

This resulting value of 6*(x) is then defined as the trough 
film thickness for the next increment. In addition, the solu- 
tion temperature distribution from the previous iteration is 
mwsed as the starting temperature distribution for the next 
increment. 

imas iterative process at each increment is repeated until 
convergence is met, and, is continued at each increment until 
the entire length of the condenser has been transversed. In- 
cremental heat rates are then summed to yield the total heat 


@eanster rate. That is: 


= * * e 
Seat al 22 EN x Q, °Ax (eqn 3.6) 
where ZFIN is the total number of axial fins. 
Once the total heat transfer rate has been determined, the 


problem is solved and pertinent data is provided as output. 


C. SMOOTH TAPERED CONDENSER SOLUTION 

The heat pipe condenser is divided into NDIV number of 
axial increments as in the finned-tapered condenser solution. 
These axial increments are then subdivided into 360 segments 


of equal length; these segments are the basic symmetric unit. 
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This unit section, is divided into a number of linear tri- 
angular elements with the same limitation as before; the 
mamber of system nodal points must not exceed 100. The 
system nodal points are then assigned x and y coordinates 
based on the input geometric parameters. 

@orstart the iterative. process, as in the finned-tapered 
case, the initial value of temperature which is an input 
Parameter is used to solve for the initial value of fin 
thickness (6*) based on the Sparrow and Gregg analysis [Ref. 
2 |. 

Once this inatial value of 6* is known, the: heat: transfer 
coefficients for the internal convective elements can be 
determined using equation (3.2). These heat transfer coef - 
£Leients are used in the Finite Element Method equations. 

The equations are solved yielding a temperature distribution. 
iesbteration 1s repeated until convergence is met, just as 
mm the finned-tapered case. 

When convergence is met, that is equation (3.3) has been 

Satisfied, a new film thickness §*(x) can be found by one of 


the following equations for §*(x): 


Bee ee Beg oneal g*)° 
Se lg See Dr, el aoe Rev, ) - 12-0 (éqn- 3.7) 
or 
gt4° 
Sh, ie] = = 0 (eqn 3.8) 
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sat 


Pet he x* cos 
u k(T, : 


Dr 


at Ty) 


Pv cos@ 


Rev 
u 


acceleration due to gravity (ft/hr7) 


0Q 
i) 


vapor velocity (ft/hr) 


< 
mM 


tae 
Ul 


local wall temperature (degrees FO 


ic 


- shear stress vapor-liquid interface (l1bf/ft*) 


Equation (5.7) defines the film thickness distribution for 
a smooth tapered rotating heat pipe derived by Daniels and 
At-Jumaily [Ref. 13]. This equation takes into account the 
Beae effects of counter-flowing vapor. Equation (3.8) is a 
modification of Equation (3.7) neglecting the drag losses. 
Equation (3.9) was developed by Ballback [Ref. 1]. This equa- 
muon also neglects drag. 

Depending on a particular control parameter which is part 
Of the input data (See Appendix B) one of these equations is 


used to solve for the film thickness (6*(x)) for the next 
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increment. In addition, the solution temperature distribution 
feeom the previous iteration is used as the starting temperature 
Gastribution for» the next increment. 

This iterative process at each increment is repeated until 
convergence is met, and is continued at each increment until 
the entire length of the condenser has been transversed, just 
moean the finned-tapered case. Total heat transfer rates are 
then determined by summing the incremental heat transfer rates 


for the entire length of the condenser by the following 


relationship: 
NDIV 
CLA tg es Q; Ax (eqn 3.10) 


D. SMOOTH CYLINDRICAL CONDENSER SOLUTION 

hoein the sSmooth-tapered case, the condenser is first 
Gavrded axially, then it is divided circumferentially into 
360 segments of equal length. These segments are the basic 
Symmetric unit section to be considered. Again the symmetric 
unit is subdivided into linear triangular elements and x and 
y coordinates are assigned to the system nodal points. 

To begin the iteration, an initial temperature estimate 
Which is an input parameter is assigned to the convective sur- 
face nodal points. Using this initial temperature estimate, 
the maximum film thickness, val When as tocated at x = 0 


max z 


ia@etermined. This maximum film thickness value is then 
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used as one of the boundary conditions in the solution of 
equation (2.18) using the Finite Element Method. Equation 


Meels)-is repeated here for reference. 


d * * 3 


8 c- I S = _2ke(Tsat - Ty)u (eqn 3.11) 


The finite element solution of Equation (3.11) provides the 
film thickness profile (6*(x)) at the midpoint of each incre- 
meme along the length of the condenser, but only the first 
increment value is applied at this point of the analysis. 

@aee S“(x) 1s known at the first increment the heat transfer 
Geefficients for the internal convective surface elements can 
be determined using equation (3.2). The steady state heat con- 
duction problem is then solved and a temperature distribution 
feeeene. unit section results. .This process is now_repeated. 

A new maximum film thickness, film profile and thus 6*(x) at 
the first increment is found based on the solution temperature 
distribution from the first iteration. With this new value for 
6*(x), the heat conduction problem is again solved for a new 
temperature distribution. This iterative process, involving 
Solution of the film profile with each iteration is continued 
until temperature convergence is met at the first increment. 
Mem convergence is met, the film profile determined on the 
iteration in which convergence was met is used to provide the 
values of 6*(x) for the remaining increments along the length 


of the condenser. 
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Using this predetermined value of 6*(x), the iterative pro- 
Sess 1S Continued at each increment until temperature conver- 
Mence is reached. As convergence is reached at each increment, 
the internal wall temperature, which is the same for all inter- 
meiewall. nodal points of the section, is stored for future 
application. 

When convergence is reached at the final increment, equation 
rid) ais once again solved for the film profile. It should be 
noted, however, that the right hand side of equation (3.11) is 
temperature dependent. It should also be noted that the temper- 
@eure varies axially along the length of the condenser. In 
Order to account for this temperature dependence and temperature 
mamtation, the right side of equation (3.11) is now determined 
for each increment using the wall temperatures that were stored 
taemedch increment. The finite element solution of the film pro- 
file will now account for the temperature variation along the 
Hength of the condenser. 

Maas final film profile provides the value of 6*(x) for each 
mierement. The iterative process of solving for heat transfer 
coefficient, temperature distribution, and temperature conver- 
gence is continued for each increment until the final increment 
is reached. 

Hquation (2.13) provides a relationship for total mass flow 
maeenas a function of position. The total mass flow rate at 


the overfall into the evaporator is given by: 
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where 6*(L) and dg*(L)/dx are the values for the film thickness 
e@a rate of change of film thickness with respect to x at the 
evaporator end of the condenser. Another relationship for mass 
flow rate based on the steady state heat conduction solution 


is piven by: 


NDIV 
ins. 560*) 2 Siee (egne3,.1 3) 
1-4 ee 


memern tact, the solution of equations (3.12) and (3.13) are 
equal, then the mass flow rate of the condensate returning to 
Eve evaporator is equal to the mass flow rate of the vapor 
Being condensed on the surface of the condenser. Or, to put 
meeanother way, continuity is satisfied. 

It should be noted that the film profile maintains the same 
Saeae shape, that is, 6*(x) at x=0 is always equal to eee and 
decreases to a specific minimum value at x=L. This being the 
aoc. at the maximum film thickness is varied, the film thick- 
Mexx profile will vary in. the same manner. For example, if 
the maximum film thickness is increased, the entire profile 
Mmiiealso increase. This will result. in an increased internal 
thermal resistance and thus lower heat transfer rate. As a 


mesult of the lower heat transfer rate, the mass flow rate as 


$5 


{atv& wes} 


determined by equation (3.13) will be less. At the same time, 
however, the greater film profile will result in a greater 
fue ofod*(x) at the overfall, x = L. Yet, since the profile 
maintained the same basic shape, the derivative at the overfall 
remains relatively constant. Thus the mass flow rate as deter- 
mined by equation (3.12) will increase. A decrease in the 
maximum film thickness will result in an opposite effect to 
the mass flow rates determined by equations (3.12) and (3.13). 
This being the case, the mass flow rates, as determined 
Myeequations (3.12) and (3.13) are now compared to determine 
eeetaeotilmyprofile is in fact the solution profile to the 
Problem. lf continuity is not satisfied, the maximum film thick- 
mess 1S varied and the entire iterative process is restarted. 
mrtoeprocessS is continued until the film profile mass flow rate, 
equation (3.12) converges towards the heat transfer mass flow 
rate, equation (3.13). When the absolute difference between 
these mass flow rates is less than a specific value, the re- 
eueting heat transfer rate is considered the solution to the 


problem. 


E. FINNED CYLINDRICAL CONDENSED SOLUTION 

As in the finned-tapered case, the condenser is divided 
axially and circumferentially into the basic symmetric section 
as shown in Figure 1. This unit is then subdivided into linear 
triangular elements and x and y coordinates are assigned to 


each nodal point. 
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Ton begin the iteratave process, as before, an initial. temp- 
erature estimate is assigned to each nodal point along the in- 
eemial convective surface. An initial trough film profile is 
determined, using equation (3.11) and this initial temperature 
estimate. In this case, however, the maximum film thickness 


(es * 


a is not.calcuiated but is. an input parameter. 


Once §6*(x) is known at the first increment, the internal 
heat transfer coefficients are determined, using equations (3.1) 
and (3.2). These values are then used in the Finite Element 
Method solution of the steady state heat conduction problem. 

A temperature distribution is determined and the iteration is 
Zepeated until temperature convergence is met. Note that a 
Mew film profile is determined for each iteration. 

As in the cylindrical-smooth condenser case, once conver- 
Bemce 15 met at the first increment, the film profile that was 
determined for the iteration prior to convergence at the first 
myerement is then used to provide the film thickness §*(x) for 
the remaining increments. 

At the final increment, a new film profile for a finned 
cylindrical condenser is then determined by solving the follow- 


ing equation developed in Chapter II: 


* 
at i es (s* e+ 6** tana) ] ae 3k CT t = Ty ue 
Oecueom hi 
£ fg 
aN eae 30 yuz* 3/4 
25% cosa | a sat avg eae 
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Pea@escription of the’ solution of this differential equation 
using the Finite Element Method is provided in Appenxix A. 
The iterative process of finding the solution temperature 
distribution for all increments is then repeated until the 
length of the condenser has been transversed. 
Mearns point, the total mass flow rate of the condenser 


returning to the evaporator is determined by the following 


relationship: 
Pew? do* 
mh = ain Fe (eot’+o*"* tana) * ZFIN (eqn: Sia5) 


where ZFIN is the number of axial fins. 
This mass flow rate is compared to the mass flow rate given 
by the following equation: 
NDIV 
pipet ZEIN*Z*. 2 ion 8a (eqn 3.16) 


i=l ae 


ust aS in the smooth-cylindrical condenser case, if the 
absolute difference between the two mass flow rates is less 
than a mass flow convergence criterion, the problem is con- 


Sidered solved. If not, § is varied and the entire itera- 


* 
max 
fame process is started again. As inthe smooth-cylindrical 
condenser, varying ae will have the same effect on the film 


profile and the heat transfer rate. Since the temperature dis- 


tribution along the condenser has been solved once and closely 
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approximates the final solution to the problem, on successive 
meerations, equation (3.14) may be used to solve for the film 
Meorile rather than equation (3.11). Equation (3.11) was used 
On the first iteration because the finite element solution con- 
verges more quickly than equation (3.14) when an estimated 
temperature is used. 

A word of caution is required. The solution of equation 
(3.14) is highly sensitive to the value of On See hf fthervini tial 
value of eee is Gineonsi stent with the actual solution, ‘eg. 
woo-smalil, the Finite Element Method solution of the film pro- 
mie will not “converge. If this is in fact the case, the prob- 
lem will be automatically terminated. A new value of ane should 
emgen be chosen and the problem restarted. 

One additional topic which should be addressed is the rec- 
tangular fin solution process. The rectangular fin profile is 
a slight modification of the finned-tapered or finned-cylindrical 
condenser solution. The only variation is that the top elements 
Or the rectangular fin are assigned heat transfer coefficients 
mame. Thus, the tip of the rectangular fin is considered 
adiabatic and no heat is transferred through this top face. 
Other than this modification, the solution techniques are the 


same as for the finned cases addressed above. 
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IV. RESULTS AND DISCUSSION 


Prior to.the heat transfer analysis of the various con- 
denser configurations, it was necessary to verify the finite 
element solution of the film profile. The development of this 
solution is discussed in Appendix A. 

Leppert and Nimmo [Refs. 8 and 9] had developed an analy- 
meedl solution, for film condensation on a horizontal plate at 
geconstant surface temperature. Their analysis and resulting 
differential equation is identical to the development in Chapter 
imeetor.a.smooth cylindrical condenser. if the acceleration due 
EO gravity is replaced by a radial acceleration term. This 
being the case, this modification was made and the analytical 
Solution was used as a reference for comparison. 

hex the test runs.of,the finite.element solution, a constant 
surface temperature was.assumed. In addition, a value for the 


) 


at the overfall were required. The value for See was determined 


maximum film thickness ona) and minimum film thickness ee 
based on a relationship developed by Leppert and Nimmo [Ref. 8]. 
The value for ba 4 was arbitrarily chosen. 

For identical geometry, surface temperature and maximum and 
minimum film thicknesses, the results of both analyses were 
identical. In order to develop confidence in the finite element 
Selution, én; was varied from 0.10 + 6% Co 0.97 =. 6 . The 


alagl max max 


Besulting profiles agreed at all locations along the length of 
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the condenser. However, when the temperature was permitted to 
mary along the length of the condenser in the finite element 
solution and a corresponding average temperature was used in the 
B@alytical solution of Leppert and Nimmo [Refs. 8 and 9], the 
profiles were no longer in agreement. This was to be expected, 
Marticularly in the case where there was a sixteen degree 
Fahrenheit variation along the length of the condenser. This 
Substantial temperature variation resulted in a significant 
waeration of fluid properties which would account for a differ- 
ence in film profile. 

Due to the agreement between the finite element solution 
and the analytical solution of Leppert and Nimmo [Refs. 8 and 
9] for a constant surface temperature, it was decided that the 
mimate element solution does provide a satisfactory representa- 
been Of the film profile. This being the case, the finite 
etement solution was then incorporated into the code to provide 
meme film profile for the cylindrical condenser. 

Once the finite element solution of the film profile was 
Werified, the heat transfer analysis could be accomplished. 

The analysis considered both copper and stainless steel con- 
densers with the following four geometries: a) tapered-smooth, 
b) tapered-axially finned, c) cylindrical-smooth and d) cyl- 
indrical-axially finned. Table I lists the geometric para- 
Meters held constant for all analyses. In all cases, the 


working fluid was water. 
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TABLE 


Condenser Geometric Parameters 
Held Constant During All Analyses 


condenser length 8.500 inches 


minimum radius OL5i575 inches 


wall thickness OQL0Si18 Iimehes 

Im addition, the following geometric parameters were also 
utilized when required. This requirement was based on the 
condenser geometry being considered, i.e., tapered-axially 


finned. 


TABLE 12 


Condenser Geometric Parameters 
Applied as Required 


0.05906 inches 


henge Or “ria 


fin nalft “angle = 20° ."56'5 degrees 
condenser cone 
half angle = 1500 degrees 


In the analysis, the heat transfer rate was determined for 
ime four different geometries listed above for both copper and 
Stainless steel. The ambient temperature was set at 60.0°F 
and the heat transfer rate was determined for each possible 


combination of the operating parameters given in Table III. 
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TABLE, i11 


Operating Parameter Matrix 


Rotational Heat Transfer Saturation 
Speed CoecfLiicrvent Temperature 
(RPM) (btu/hr-ft?-F) (degree F) 
700.0 OO (0) 90.0 

1400.0 500.0 Z0'20 
2800.0 1000-0 ISO 0 
ESO 20 


Thus, for each condenser geometry, there was a total of 72 


analyses, 36 


for the case of the condenser with a copper wall 


and 36 for the case of the condenser with a stainless steel 


wall. 


The first condenser geometry considered was a smooth con- 


memser. Figures 8-13 compare the heat transfer rates of smooth 


cylindrical condensers with those of smooth tapered condenser. 


mieparticular, Figures 8,9, and 10 indicate the results of the 


ematyses of smooth copper condensers at rotational speeds of 


700, 1400 and 2800 revolutions per minute(RPM) respectively. 


Figures 11, 12, and 13 are for smooth stainless steel conden- 


mems at 7/00, 


1400, and 2800 RPM respectively. For both stain- 


less steel and copper smooth condensers, the following general 


observations 
Eeciticient, 
condenser is 


the external 


apply: a) For the same external heat transfer 
the heat transfer rate for the cylindrical smooth 
less than the equivalent tapered condenser. bb) As 


heat transfer coefficient increases, this difference 
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sauheat transfer rate also,increases.. .c). The rotational speed 
Waseda greater effect on the tapered condenser, heat transfer 
rate. For example, the maximum heat transfer rate of a copper 
tapered condenser will increase by a factor of 1.67 when the 
mocational speed is»increased from 700 to 2800 RPM. In. the 
Szlindrical copper condenser, for the same change in rotational 
Speed, the heat transfer rate only increases by a factor of 1.51. 

These same observations hold true for the smooth stainless 
Steel condensers. But, due to a greater thermal resistance in 
the wall of the stainless steel condenser, the heat transfer 
matees-are less forall cases considered. It should be noted 
that the thermal conductivity of stainless steel is only 4% the 
eeermal conductivity-of copper. This accounts for the increased 
enermal resistance. 

Bigunes (14;-15, and*:l6 compare axially, finned cylindrical 
With axially finned tapered copper condensers at 700, 1400 and 
BeeOeRPM:- respectively. -Note that the heat transfer rates of 
Bae cylindrical condensers are only slightly less than those of 
Pee tapered condensers. This is because the heat is primarily 
meemorerred through the extended surface, »i.e., the fin. Thus 
the film condensate in the trough has less effect on the heat 
Peanster rate. In fact; the average difference in heat transfer 
rate for 700 RPM and an external heat transfer coefficient of 
mo@ebtu/hr-ft*-F is 12.65%; This difference increases to 15.6% 
as the heat transfer coefficient is increased to 1000 Btu/hr- 


ft*-F. However, as the rotational speed is increased to 2800 
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RPM, the corresponding average differences in heat transfer 
rates decrease to 12.58% and 13.33% respectively. 

Fapuves 175008" and*19 compare the heat transfer rates of 
axially finned stailness steel cylindrical condensers with 
faese Of axially finned stainless steel tapered condensers 
mee od, 1400 and 2800 RPM respectively. Note the difference 
ee neat transfer rate increases as the external heat transfer 
Beerficient increases more than in the case of the copper 
Condensers. At the low heat transfer coefficient, the limit- 
Ms 


However, as the external heat transfer coefficient is increased, 


me thermal resistance is that of the external surface ple bepets 
the limiting thermal resistance becomes that of the wall due 
forthe low thermal conductivity of stainless steel. 

Another) observation to be noted is the fact that the rota- 
mremal) Speed has a greater effect on the heat transfer rate of 
the cylindrical condensers than on the tapered condensers. As 
Ene rotational speed increases, the film thickness decreases 
due to the greater centrifugal force exerted on the film. For 
the cylindrical condenser, the maximum film thickness is much 
Eereater than for a tapered condenser. For example, an axially 
finned stainless steel cylindrical condenser rotating at 1400 
RPM with an external heat transfer coefficient of 1000 Btu/hr- 
ft?-F has a maximum film thickness twice that of a corresponding 
eapered condenser. This being the case, higher rotational 
Speeds will have a greater effect on the greater film thickness 


ema the difference in heat transfer rates will decrease. 
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Pegures 20), 21, and 22; compare the heat transfer rates of 
copper smooth cylindrical condensers with copper axially finned 
cylindrical condensers at 700, 1400 and 2800 RPM respectively. 
Time figures indicate that for low external heat transfer coef- 
Mecwents,, little:-is to be gained by the addition of axial fins. 
But, as the external heat transfer coefficient increases, the 
advantage becomes significant. As an example, consider Figure 
21. For h=100 Btu/hr-ft?-F, axial finning increase the heat 
m@eanster rate by 16%. But, for an external heat transfer coef- 
ficient of 1000 Btu/hr-ft?-F, the heat transfer rate increases 
by 194%. 

Note also, that as rotational speed increases, the advan- 
tage to be gained by axial finning decreases. This canbe 
Smprained by the fact that for an axially finned condenser, 
the majority of the heat is transferred by the fin surface. 

As the rotational speed increases, the film thickness in the 
meeugh decreases. This+in turn will expose slightly more fin 
Surface area. On the other hand, for the smooth condenser, 
the decrease in the film thickness will have a greater effect 
on heat transfer rate in that the thermal resistance of the 
film has decreased. Comparing the two geometries, the change 
in overall thermal resistances for the smooth condenser will 
be greater than that for the axially finned condenser account- 
ing for the slight decrease in advantage with increasing 


rotational speed. 
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Breures 25, 24; and Z5 correspond to Figures 20, 21, and 
mepuc tor stainless steel condensers. The results are similar 
fo the copper situation discussed above, by the effect caused 
mye the increasing external heat transfer coefficient is not 
as dominant due to the high thermal resistance of the stainless 
steel wall material. 

Figure 26 compares a smooth cylindrical copper condenser 
with a smooth cylindrical stainless steel condenser at 1400 
mm. AS to be expected, the heat transfer rate of the copper 
Bemdenser 15 greater than that of the stainless steel condenser 
due to the difference in the thermal conductivity of the two 
Materials. The difference in the heat. transfer rate is least 
fon 2a low external heat transfer coefficient where the external 
eeermal resistance is dominant. As the heat transfer coeffi- 
cient increases, the thermal resistance of the wall of the 
Cendenser becomes more important resulting in an increasing 
difference between the two condensers. 

Figure 27 provides a comparison of axially finned cylindrical 
copper and stainless steel condensers at 1400 RPM. The advantage 
of copper over stainless steel is obvious. Note that the heat 
transfer rate for the copper condenser at 500 Btu/hr-ft?-F is 
nearly identical to the stainless steel condenser heat transfer 
mace at 1000 Btu/hr/-ft2-F indicating the advantage of copper 
Over stainless steel. 

The final analysis compared the heat transfer rate of an 


axially finned copper condenser with a triangular fin profile 
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forthe heat’ transfer rate of an axially finned condenser with 

a rectangular fin profile. The rectangular fin was assumed to 
have an adiabatic tip. This comparison was accomplished for 
both cylindrical and tapered condensers. Table IV lists the 
parameters used in the analysis. These parameters are in addi- 
Ewon to those listed in Table I. This comparison was only 
accomplished for the one set of operating parameters listed in 
Mabie IV. Note that the operating parameters chosen were the 


median values. 


TABLE IV 


List of Parameters Used in 
Rectangular/Triangular Fin Profile Comparison 


Heat Transfer Coefficient SOs hrort ose 


Rotational Speed 1400 RPM 


Saturation Temperature 120.0 degrees F 
Bin Height = 0.05906 inches 


Mable V lists the results of the analyses. 


TABLE V 


Results cf Rectangular/lrianeular 
Ein Profile Comparison 


CONDENSER GEOMETRY FIN PROFILE Vatetu/hr) 


tapered triangle 6168493 
rectangle G1:43 31013 
cy hindrircal triangle 5298:.4 
rectangle 5338.00 
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Note that the heat transfer rate varied by only 0.4% for both 
fae tapered and cylindrical condensers. The reason for this 
misteniticant variation is provided in Table VI and VII. 

Table VI lists the convective surface temperature distribution 
of a symmetric unit section of a tapered condenser at the middle 
increment of the condenser. Temperature location 1 is located 
Gemeeneitip of the fin. For the rectangular profile, temperature 
Heedation 1 is located at the intersection of the adiabatic 
pumrace (the tip of the fin), and the vertical surface of the 
mine Lemperature location 5 is located at the base of the fin. 
Locations 6-9 are located in the trough region and the remaining 
memperature locations are along the external surface of the 
weertion. Thus temperatures 1-5 provide the temperature, distri- 
bution along the convective surface of the fin. Note that the 
tPemperature distribution is lower for the rectangular fin. 

mags, ithe driving force for heat transfer, the temperature 
difference between the saturation temperature and the surface 
Mmemmenmature 15 greater for the rectangular fin. ‘Thus, in spite 
Sepene fact that the fin surface area for heat transfer has de- 
creased by 11%, the average surface temperature difference has 
increased by 25%. It should also be noted that the average 

heat transfer coefficient for the fin surface for the rectan- 
gular fin also decreases by 12%. However, the increase in temper- 
ature difference is the dominant change that only allows a 0.4% 


@ecrease in heat transfer rate. 


69 


od 


i i & 


——- va ww 
dzot yp? th vine Ye bee 3 
ies p. 
e277 You BOP ROT at -? 7 SeE5 Bros y 
pre I! Oidel i Basevene rae 


Nd 
ms3 232 tHe SVERs 
'o noLtgosers 


1s2nsbnag 


TABLE VI 


Surface Temperature Distribution for Rectangular 
and Triangular Axially Finned Tapered Copper 
Condensers at Middle Increment of Condensers 


Temperature Triangular Fin Rectangular Fin 
Location Temperature (F) Temperature (F) 
1 TS, 58S 7, 5 
2 118.45 17, ASS 
3 118.04 R728 
4 17, 160 TG, $36 
5 TAG? OS TRG, 259 
6 T16..282 D166. 2455 
7 P16, {Gi IEG: 25'S 
8 TRG, 2519 D6, £35 
9 116.48 BG. 354 
10 16, S47 TG, 255 
aa 116.46 ING... F2y5 
12 116.46 TSG, 7231: 
1S 116.44 TRG, 3210 
14 IEG, £452 IG 157 
155 D6 350. APISG, 5122 
16 pLaWione 25s i 6, 007 
7 DNG, 2277 116.04 
18 TG: 2216 156, 2033 
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TABLE VII 


Surface Temperature Distribution for Rectangular 
and Triangular Axially Finned Cylindrical Copper 
Condensers at Middle Increment of Condensers 


Temperature Triangular Fin Rectangular Fin 
Location Temperature (F) Temperature (F) 
id 48 259 Ae? 25555 
2 1809 117.44 
5 62 ils areas 
4 17 SO eGR Se7 
5 £16.55 116.36 
6 11.6 558 LPOwcs 
ul TO; 208 £16.07 
8 11 6).22 LEG... 1 
9 116-20 Laon 9 
10 LAG 0:3 115.80 
ite 1160/2 1S ...79 
12 115.07 TAS eyo 
a3 1ESR 99 ELS 77 
14 115.97 LS e7S 
15 5.94 Gi Een 
16 IDES el LS 267 
17 115.88 145-20 7 
18 115. 88 115 .66 
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Taple Vil indicates a Similar Situation exists in the 
Emlindrical condenser. Again note the decrease in the average 
Surface temperature of the fin for the rectangular profile. 
The surface area and heat transfer coefficient has decreased, 
but the increased temperature difference compensates for these 


Changes, limiting the overall decrease in heat transfer rate. 
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Figure 8. Heat Transfer Rate versus Saturation Temperature 
fOr somooth Cylindrical and Tapered Copper Con- 
densers at 700 RPM 
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Figure 9. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Copper Con- 
densers at 1400 RPM 
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Figure 10. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Copper Con- 
densers at 2800 RPM 
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Figure 11. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Stainless 
steel Condensers at 700 RPM 
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Figure 12. Heat Transfer Rate versus Saturation Temperature 
: for Smooth Cylindrical and Tapered Stainless 
Steel Condensers at 1400 RPM 
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Heat Transfer Rate versus Saturation Temperature 
for Smooth Tapered and Cylindrical Stainless 
Steel Condensers at 2800 RPM 
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Figure 14. Heat Transfer Rate versus Saturation Temperature 
for Axvaliy Finned Cylindrical and Tapered Copper 
Condensers at 700 RPM 
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Figure 15. Heat Transfer Rate versus Saturation Temperature 
for Axtally. Fimned .Cylandrical..and. Tapered, Copper 
Condensers at 1400 RPM 
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Figure 16. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Copper 
Condensers at 2800 RPM 
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Figure 17. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Stain- 
Mess steet Condensers at 700° RPM 
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Figure 18. Heat Transfer Rate versus Saturation Temperature 
fom Axially Fanned Cylindrical and Tapered Stain- 
less Steel Condensers at 1400 RPM 
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Heat Transfer Rate versus Saturation Temperature 
for Axially Finned,. Cylindrical and Tapered Stain- 
less mreele Condensers! at 2800 RPM 
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Beeore 20. Heat Transfer Rate versus Saturation Temperature 
forssmooth and, Axiaillly Finned Copper Cylindrical 
Condensers at 700 RPM 
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Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Copper Cylindrical 
Condensers at 1400 RPM 
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Figure 22. Heat Transfer Rate versus Saturation Temperature 


for Smooth and Axially Finned Copper Cylindrical 
Condensers at 2800 RPM 
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Figure 23. Heat Transfer Rate versus Saturation Temperature 
for smooth and Axially Finned Stainless Steel 
Gydandrical ‘Condensers at 700’ RPM 
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Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Stainless Steel 
Cylindrical Condensers at 1400 RPM 
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Figure 25. Heat Transfer Rate versus Saturation Temperature 
fOr usmeoth and Axially “Finned ‘Stainless’ Steel 
Cylindrical Condensers at 2800 RPM 
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Preure 26. Heat Transfer Rate versus Saturation Temperature 
for Copper and Stainless Steel Smooth Cylindrical 
Condensers at 1400 RPM 
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Figure 27. Heat Transfer Rate versus Saturation Temperature 
for Copper and Stainless Steel Axially Finned 
Cylindrical Condensers at 1400 RPM 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

the heat transier rate of a.cylindrical condenser is less 
than an equivalent tapered condenser. This decrease in heat 
@eanster rate is most significant in a smooth condenser, where, 
depending on the external heat transfer coefficient and rota- 
tional speed, can be as great as 45%. This decrease in heat 
Granstfer rate becomes less significant for an axially finned 
Gemcenser where the average decrease, for the range of heat 
muansrer coefficients examined, was 15% for a copper condenser 
gastos, £0r a stainless steel condenser. When such factors 
mmeeme cost and difficulty in manufacturing tapered axially 
finned condensers are considered, this 13% decrease in heat 
feanster rate becomes tolerable. From a practical standpoint, 
development and analysis of cylindrical axially finned con- 
densers should be encouraged by these results. If environ- 
mental conditions permit, copper should be preferred over 
Geainiess steel due to its exceptionally high thermal con- 


ductivity and resulting higher heat transfer rate. 


B. RECOMMENDATIONS 
I) Build and experimentally test both smooth and axially 
finned cylindrical condensers to obtain experimental data 


for comparison with results of this analysis. 
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2) Develop models for rectangular and trapezoidal fin 
profiles with non-adiabatic tips and incorporate into code. 

3) Experimentally test axially finned cylindrical con- 
densers with rectangular and trapezoidal fin profiles to 


obtain data for comparison with theoretical results. 
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APPENDIX A 


FILM PROFILE. FINITE ELEMENT SOLUTION 


A. SMOOTH CONDENSER 
mae analysis of the film profile in a smooth cylindrical 
condenser developed in Chapter II resulted in the following 


Srainary, nonlinear, second order, differential equation: 


d dg* ate SkelTeaemen ty) u 

* * a a W 

oa. ia (OC! igcch ages (eqn A.1) 
Pp wr he, 


firs equation can be rearranged and expanded to yield: 


2 


AS d2 <= dg* 


6 + 5675 ( ) = -K (eqn A.2) 
dx? dx 
where K = os Gees - ieee 
eee. 7 
Qe oe 


The statement of the problem for the formulation of the 


Finite Element Method is: 


Buedog* 43 a are S 
6 eee) + 36 (Gee) K (eqn A.3) 


meet x = 0, g* = cee 


="05 da */dx = 0 
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The domain (length of the condenser) is divided into ele- 
ments Of length Ax with the exception of. the, first, and final 
elements which have length Ax/2, where Ax is the length of the 
condenser divided by the number of axial increments (NDIV). 
mesystem nodal point is located at each end of an element. 
Tus, system nodal point 1 is located at x=0 and the last sys- 
tem nodal point, which is equal to the number of elements plus 
irs located at the overfall into the evaporator, x=L. All 
iaternal nodal points are located at a position corresponding 
mene Midpoint of the axial increment. 

Herine. the approximate value of the film thickness in the 


following manner: 


h 

* 

6 2°, Geis -oUG.d-.=)G d (eqn A.4) 
ah 


il 
ct 
» 
oO 


where 6 approximate value of the film thickness (6*). 
oe the global basis functions... 
n = the number of system nodal points. 


d,=-the solution; vector. 


On an element level, equation (A.4) becomes: 
e ue 4 
be = 0 (x) = » God. (eqn A.5) 


where cele the locals basis> functions. 
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Define two degrees of freedom at each nodal point, i.e., 
both 6 and dé/dx are continuous. Thus, the local basis func- 


tions used in the finite element solution are: 


eprn hee st: ¢ nS (eqn A.6) 

- ae = ' (eqn A.7) 

g, = B - B (eqn A.8) 

Se = BS i (eqn A.9) 
where & = the length of an element. 


Zi and Bee magnitude basis functions. 
g, and g, = slope basis functions. 
Substituting the approximate film thickness OF into the 
@yviterential equation (A.3) results in: 
4 2x 223 2 
mote fe. 35 “(t8ey = -x (eqn A.10) 
e : e 
dx dx 
To remove the "nonlinearity" from the problem, equation (A.10) 
is modified inthe following manner: 
dé 


= + 3n?n! —== -K, (eqn A.11) 
dx? dx 
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n is defined as the approximate value of the film thickness 
arom the previous iteration. In like manner, n' is defined 
as the approximate value of the rate of change of film thick- 
fessowith respect to x from the previous iteration. 


rerming the residual of equation (A.11) yields: 


Ra - aes dja aS COeat aK (eqn A.12) 


Invoking the Galerkin criterion for the determination of the 


Solution vector d, i.e. 


fG;R,dx TOTES SI A yee Sc | (eqn A.13) 
yields: 
ere ah a ie z 
meecG ddx + 3n-nesiGG ddx + KfGdx = 0 (eqn A.14) 


Each of the integrals in equation (A.14) are defined in the 


following manner: 


" 


T n 
[GG ddx = |l™ eg 
~ ~ lh 


di dx (eqn A.15) 
mee didx =| / j\€. y didx (eqn A.16) 


This integration results in the following 4x4 local elemental 


A matrix for any element: 
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ima similar manner, let 


ce 


n ! 
HGG ddx = & fgg d dx (eqn A.17) 
l 2 2 


and on an elemental level: 


Be <87 g gz a > 
T 32 
fogg d dx aig a 83 d dx (eqn A.18) 
2 g) 


This integration results in the following 4x4 local elemental 


Bemiatrix for any element: 


-1 Q 1 -2 
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Pastly, Let 
n 
Gdx- =*"E"*s\ gdx (eqn -A.19) 
a 
Or, on an elemental level this becomes: 
fogdx = fy g3 dx (eqn A.20) 


This integration results in the following local elemental 


column F vector: 


(FI, = 


muus, £0r a given element, equation (A.12) becomes: 


MeeAl< to3in5° im,’ (Ble) dg = -X, FZ (eqn A.21) 
As mentioned above, n and n' are the approximate values of §* 
and dg*/dx respectively from the previous iteration. For the 
m@@etial iteration; n is set equal to bak andann iS \yset equal 
Bo 0. 

Each elemental matrix is placed in a global system [A] 


Matrix with the location in the matrix based upon the local/global 
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modal point correspondence. For example, for element Nr. 2 

fhe nodal points corresponding to local nodal points 1,2,3, 

and 4 are global nodal points 3,4,5, and 6 respectively. 
Therefore, the sum of the two 4x4 elemental matrices ([A] and 
[B]) and multiplying constants will form a single 4x4 local 
maeecin,. 1. Element (1,1) of this local matrix will be placed in 
the global A matrix location (3,3) and element (4,4) of the local 
elemental matrix will be placed in the global A matrix location 


(6,6). Ultimately, the following global system would be assembled: 


[A] . d = F (eqn A.22) 


Note that the global A matrix would be an mxm size matrix 
Where m is equal to twice the number of system nodal points to 
meount,+0r the two degrees of freedom,at each nodal point. 
Beaecimilar fashion, d and F would be column vectors of size m. 

Once the system is assembled the boundary conditions are 
@pplied. This is done in. the following manner: For boundary 


eendition (a), A(1,1) is set equal to 1.0 and the remaining 


Paements in the first row, i.e. A (1,1).1=2,3,4,...m are set 
i k % Vie ieee 
eae, tO 0.0. Then F(1) is set equal.to jas: Snax iS imi 


tially determined by a relationship developed by Leppert and 
Nimmo [Refs. 8 and 9]. This establishes the value of d(l1) as 
[a oil assamilar manner, boundary. condition. (b),is applied 


by setting A(2,2) equal to 1.0 and all other second row elements 


G@eethe global A matrix are set equal,to 0.0. Then, F(2) is set 
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Pamala to 0.0), . Thais establishes the solution vector d(2). as 0.0. 
moce that d{2) corresponds to the slope at the first nodal point. 
One additional boundary condition is required; this is the 
Matue of the film thickness at x=L. This boundary condition is 
Mecessary to completely specify the problem. The value of the 
film thickness at the overfall may take on any value depending 
emethe geometry at the overfall. In the case of this analysis, 


Pans value was taken as 0.25 + 6§* this value was chosen for 


max” 
the following reason: Leppert and Nimmo [Refs. 8 and 9], in 
a Similar analysis for laminar film condensation on a horizontal 
Surface derived an analytical solution to equation (A.1), assum- 
me a constant surface temperature. They found the film profiles 
for 6* with the overfall value less than 0.40 - ae os Were \eSsen- 
tially constant and thus any value of §* at the overfall less 
than 0.40 - Sten would result—in the’ same profile. In the veri- 
fication of the finite element solution, not only was this found 
bo be the case, but it was also found that the heat transfer 
meee Was relatively insensitive to the shape of the film profile. 
Mieract, it was found that the film thickness at the overfall 
could be increased to a value as great a 0.90. oe and the 
meoulting variation in heat transfer rate was only\10%. This 
being the case, the value of 0.25 - See was arbitrarily chosen 
fOr an” 
The third boundary condition was applied by setting A(m-1,i)= 
0.0 where i=1,2,3...m. Then, the global matrix element A(m-1,m-1) 
was set equal to 1.0. Finally, F(m-1) was set equal to 0.25°5 75 
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Once ali three boundary conditions were applied, the sys- 
tem given by equation (A.22) was solved for d. The values of 
Pee inpcoximatce cLidmathickness., (|tse2, d(1), f=1,55455...m-1 are 
then compared to the values of film thickness from the previous 
meeration. (Ifthe relative difference is less than or equal 
Gena specified convergence criterion (i.e., 0.0001) at all 
meaad points, convergence 15 met and the latest values of d 
are the solution values of 6*. 

tee convergence 1S not met, the values of d are saved for 
tne next iteration where they are. used to determine n and n' 
as discussed above. This iterative process is continued until 
convergence is met or until a maximum number of iterations have 


occurred. 


B. AXIALLY FINNED CONDENSER 

The finite element solution for the film profile in an 
a@xdally finned cylindrical condenser is very similar to that 
Saeea Smooth cylindrical condenser. For this reason, only the 
variations in the development will be addressed. From Chapter 
Meeecae differential equation for mass flow rate in an axially 
finned cylindrical condenser is given by: 
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This equation can be rearranged and expanded to yield: 


d76* Pan 5 dg* ea dor x4 dg * 
aa (ec#"s0.6*"° tana) + + (Seo 7 ao tana a 


a: is * 
= Ky K,6 (eqn A.24) 
where R= eel eat TYJHE 
Pe wr Reg 
fi 
Ak «(T To )uz* | 
K,= 2cosa ~ oo ove 


Substituting equation (A.5) into equation (A.24) results in: 


ae. _thickgess, dé e006, 4 ds. 
oF (ya 6, tana) + (3e6. SF we 45, tana =) 
mee Kae K,5, (eqn A.25) 


Hoeremove the "nonlinearity" from the problem, equation (A.25) 


is modified in the following manner" 


BCS. dd, 
(ey*+y>°tana) + 
dx dx 


(Sey y'+ 4y“*y"' tana) 


me Fe Koh (eqn A.26) 
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Here, y and y' are defined by the following relationships: 


y = 85 + R*(65 Oa) (eqn .Al2:7) 
yi= S, + R*(G, “G,) (eqn A.28) 
where 8, ="approximate value of the film thickness for 
the present iteration. 
Sobel approximate value of film thickness from the 
previous iteration. 
Sc = approximate value of the derivative of the film 
thickness for the present iteration 
ci approximate value of the derivative of the film 
thickness fromthe previous iteration,. 
Re = relaxation factor. 


imese two variables, y and y', are in actuality, adjusted approx- 
aiation of film thickness and derivative of film thickness re- 
Seeetively. This adjustment is required in order to converge to 
a solution. 

The finite element solution is now identical to that of the 
Smooth condenser, that is, the residual is formed, the Galerkin 
criterion is invoked, and identical 4x4 local matrices are 


derived. Finally, the equivalent of equation (A.21) is formed. 


4 5 Sarat ut 5 
Mey. + Y, tan) PN Ge uy Y. tane)i[B | "*d. 


(eqn A.29) 
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Notice that this equation has two forcing terms. The additional 
term (K5Y) resulted from the nonlinearization of the problem. 
must an in the smooth condenser film profile solution, the 
global system given by equation (A.22) is formed, the boundary 
Eougitions» applied and’ the system solved for a solution vector 
Memebiterative process is continued until convergence is met. 
With each iteration, y and y' are updated and used for the 
next iteration. When convergence is met, the latest values of 


mie i=1,3,5,...m-1 are the solution values of 6*. 
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APPENDIX B 


USER'S MANUAL 


this-appendix-describes the data-cards required to use 
the computer code. 

iimcedata.is.divyided-into. “biocks"™.for-convenience:~» Each 
@aee of this user's guide is a separate block. For each block, 
a general description, the required format, and appropriate 
comments are provided. 

ie is imperative that input data be consistent or errors 
will result. For example, if a smooth geometry is being analyzed, 
the finite element parameters must also result in a smooth model. 
imeaddition, all data fields must be filled with an input value, 
even 1f that value 1s not needed for the analysis. For example, 
iid smooth analysis, no fin half angle is required for the cal- 
Sutations; however, a value of the correct format must be pro- 
muaed in the fin half angle field or an INPUT/OUTPUT error will 


result. 
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DATA BLOCK A 


DESCRIPTION: FINITE ELEMENT PARAMETERS 


FORMAT : 


8I5 


NCMREC NCMTRI NRWFIN NRWTRF NCMTRF NCOL NPRNT 


CONTENTS 


NDIV---Number of axial increments. Must be less than 
Or equal £o 50. 


NCMREC-Numbex Of, columns of finite.eclements in. the rec- 
Eanoular section of fin. May be equal to 0 if 
fen 1s, triangular only. 


NCMTRI-Number of columns of finite elements in tri- 
angular Section of tin. May be equal to 0 if 
only rectangular fin. 


NRWFIN-Number of rows of finite elements in the fin 
section of model. May be equal to 0 if a smooth 
condenser. 

Note: If triangular or trapezoidal fin, NRWFIN 
must equal NCMTRI. 


NRWTRF-Number of rows of finite elements in wall section 
of the model. 


NCMTRF-Number of columns of finite elements in the 
trough section of finned model. Set equal to 
0 if a smooth condenser. 


NCOL---Total number of columns of finite elements. 
Must be equal to NCMREC+NCMTRI+NCMTRF for a 
finned model. 


NPRNI--Print control number. 
[If equal to 1 correspondence table and major 
elements of finite element model will be printed 
out . 
ie equal to.0; output will be suppressed. 
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DATA BLOCK B 
DESCRIPTION: FLUID, FIN MATERIAL SELECTION PARAMETERS 


FORMAT: 215 


—|—e—-]|— eg ge Me eg ge ewe Be MX Me ee ewe ee eM ee eB Ke ee ewe ee ee ee Ke ewe ee ew ee Be Ke ewe ee ee ee ew ee ee Ke eK ee eee eK 


FIELD CONTENTS 


i fFEULD-Tfméquad toe, working*fluid is water. 
if equal to 1; working. fiuad is freon. 

Z IFIN---If equal to 0, condenser wall material is 
copper. 


If equal to 1; condenser wall material is 
Sieaimless steed. 
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DATA BLOCK C 
DESCRIPTION: CONDENSER GEOMETRY 


FORMAT: 6G10.5 


CLI REASEL THICKI BFINI CANGL FNWTHI 


FIELD CONTENTS 


d! CLI----Condenser length (inches). 


2 RBASEI-Inside radius to wall of condenser at condenser 
end (inches). 


3 THICKI-Wall thickness (inches). 


4 BFINI--Fin height (inches). Must be set equal to 
0.0 if smooth condenser. 


5 CANGL--Condenser half angle for tapered condenser 
(degrees). Must be set equal to 0.0 for 
cylindrical condenser. 


6 PNWIHL-Width of rectangular portion of trapezoidal 


Or rectanpular fin. (inches). Lf triangular 
tin, set equal to 0.0. 
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DATA BLOCK D 
DESCRIPTION: INTERNAL FIN GEOMETRY 


FORMAT: 2G10.5 


FIELD CONTENTS 


if FANGL--Fin half angle (degrees). Set equal to 0.0 
for*’smooth condenser or rectangular fin. 


2 ETOEO--Ratao or trough width to fin base width. 
Determines spacing between’ fins. 


je Tail 


DATA BLOCK E 
DESCRIPTION: CONVERGENCE CRITERIA 
PORMAT: 2G10.5 


FLELD CONTENTS 


1 GiGl=——Tempegarare sconviergence icriterien. ™ Usied ito 
determine solution of two dimensional steady 
state heat conduction problem. 


2 GRITDL-Mass flow convergence criterion. Used only 
in cylindrical condenser analysis for mass 
flow convergence test. 
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DATA BLOCK p 
DESCRIPTION: OPERATING PARAMETERS 


FORMAT: 5G10.5 


RPM HINF TINTL TSAT TINF 


PIELD CONTENTS 


1 RPM----Rotational speed (revolutions per minute). 


Z, HINF---External heat transfer coefficient 
(Btu/hr-ft -deg F). 


5 TINTL--Initial temperature estimate (degrees F). 
4 TSAT---Saturation temperature (degrees F). 
5 TINF---Ambient temperature (degrees F). 


eS 
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DATA BLOCK G 
DESCRIPTION: OUTPUT PRINT CONTROL 
FORMAT: 415 


IUNITS NFLAG1 NFLAG2 NFLAG3 
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FIELD CONTENTS 


1 IUNITS-Output units control number. 
Ri LUNLTS. = 0... .caleulated results will be pro- 
vided.in English units. 
ii LUNTTS = 1. anput parameters will be repeated 
im ot “Unies, and icalculaced results will be pro- 
Wided, din ol units. 
ff [UNITS = 2, input parameters will be repeated 
im ol, units and output rTesult’s will be provided 
ijt) Deena een? bush, and Sil units . 


2 NFLAG1-The first axial increment at which the parameters 
listed, under remarks will be provided as output. 
3 NFLAGZ2-The final increment at which the parameters 
listed under remarks will be provided. 
4 NFLAG3-The step change in increments between NFLAG1. 
REMARKS 


1) No matter what value of IUNITS is used, input parameters 
felt aiways first appéar in English units. 

2) For increments indicated by NFLAG values, the following 
parameters will appear: a) x-coordinate, y-coordinate and 
temperature at each nodal point, b) length of element, heat 
transfer coefficient and heat rate per unit length for each 
convective boundary element. 

a) AS 2 minimum, the values of 1, 1, 1 must be provided as 
input for NFLAG1, NFLAG2, and NFLAG3 respectively. 
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DATA BLOCK H 
DESCRIPTION: TAPERED CONDENSER SOLUTION METHOD 


FORMAT: 115 


FIELD CONTENTS 


1 NSOLVE-Tapered solution control number. 
Form tapered, axially finned condenser. ‘set 
NSOLVE=1. 


For tapered smooth condenser, NSOLVE must be 

set to one of the following three values: 

Seto NoOlLVE = 2 if solution of film thickness is 
to be based on Ballback's {[Ref.° 1] 
equation. 

set, NSOLVE. = 3 if solution of film thickness is 
to be based on Daniels and Al-Jumaily [Ref. 13] 
equation, neglecting drag terms. 

Set NoOLVEL— 4 1: “Solution of. film? thickness: 15 
to be based on Daniel's and Al-Jumaily [Ret. 
13] equation, with drag effects included. 


REMARKS 


1) NSOLVE is only used in tapered condenser analysis. 
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DATA BLOCK I 
DESCRIPTION: CYLINDRICAL CONDENSER ANALYSIS PARAMETERS 


PORMAT: 315,2G615.10 


—e ee ee ee ee ew er ewe ee ee ee eee ew ee ew ewe ew ew eee ew Kee ee ew eM ew eK ewe eM Ke eM er eer er er er er er er er reer eee 


NONCE ITERMX. ITRPRT RELAX DELMAX 
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FIELD CONTENTS 


1 NONCE=-=Sangle iteration parameter. 
i: NONCEW= 1, oniysomeviteration will, be 
permitted. 


i NONCEs= 0, .2terataons wall continue until 
convergence or maximum number of iterations 
is reached. 


2 ITERMX-Maximum number of iterations permitted in analysis. 
3 PIRPRioLeeration print @control ‘parameter . 
iL. lTReRL, = i. mass, slow convergence test results 


will be provided for each iteration. 
ff—I1TRPRE = 05 masisi.flow econvergence test results 
Will only be provided on final iteration. 


4 RELAX--Relaxation variable used in finite element solution 
od Gylandrical finnedwfilm profile. 


5 DELMAX-Initial estimate of maximum film thickness used in 
SolUGBOn, Of, Cy lindrtcal finned »film,profile: 


REMARKS 

1) Above parameters are only used in cylindrical condenser 
analysis. 

2) Recommended value of RELAX is 0.80. It is sometimes necessary 


Pomadjust this value plus or minus 0.05°to reach film profile 
convergence at small film thickness values. 

Bye input value of DELMAX is only used in cylindrical finned 
analysis. For the cylindrical smooth condenser, DELMAX is 
internally generated. 
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APPENDIX C 
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